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SUPERSONIC  COMBUSTOR  INSULATION 
ABLATION  ANALYSIS  AND  TESTS 

Robb  W.  Newman  and  Harold  G.  Fox 
Johns  Hopkins/Applied  Physics  Laboratory,  Laurel,  Maryland 

ABSTRACT  /  j/  H 

The  design  of  passively  cooled,  hypersonic  tactical  missiles  poses  many 
severe  structural  problems  whose  solutions  are  beyond  the  current  state-of-the- 
art.  Specifically,  the  heating  conditions  in  the  supersonic  combustor  of  such 
a  missile  approach  those  in  the  throat  section  of  rocket  motors  but  are  more 
damaging  to  most  materials  because  the  flow  is  oxidizing  rather  than  reducing. 

Carbon-carbon  materials  are  being  considered  for  the  supersonic  combustor 
wall  structure  because  they  have  attractive  strength-weight  characteristics  at 
elevated  temperatures.  Since  these  materials  have  poor  oxidation  resistance, 
metal  loaded  carbon-carbon  materials  which  form  an  oxidation  resistant  layer  are 
also  being  considered.  In  order  to  design  such  a  combustor,  the  designer  must 
be  able  to  accurately  predict  the  erosion  rate  of  candidate  materials  so  that  a 
realistic  balance  can  be  achieved  between  weight  and  performance.  This  paper 
addresses  the  problem  of  thermochemical  erosion  of  graphitic  materials,  the 
development  of  procedures  for  the  analytic  prediction  of  these  erosion  rates 
in  a  combustor  environment,  and  experin^ntal  procedures  used  to  validate  the 
analytical  model  and  evaluate  candidate  combustor  liner  materials. 


INTRODUCTION 

Previous  reviews  and  evaluations  of  graphite  ablation  theory  and  experimental 
data^»^  have  been  conducted  for  blunt  body  atmospheric  reentry  analyses  as  part  of 
the  Aerospace  Nuclear  Safety  Program.  Although  these  references  are  slanted  to  en¬ 
vironmental  conditions  and  configurations  entirely  different  than  those  under 
consideration  herein,  the  basic  theory  and  simplified  model  are  applicable.  There¬ 
fore,  it  was  the  objective  of  this  study  to  modify  the  existing  model  in  order  to 
consider  the  oxidation  of  graphite  wall  materials  in  the  presence  of  a  supersonic 
combustor  environment  during  a  flight  profile. 

During  the  powered  flight  portion  of  the  ramjet  trajectory,  the  fuel  is  burned 
in  varying  mass  ratios  with  air,  releasing  heat  and  by-products  in  the  combustor 
chaiid}er.  The  resulting  high  temperature  ^as  mixture  flows  through  this  plpe-llke 
configuration  at  supersonic  velocities  while  transferring  heat  to  the  combustor 
wall.  The  heat  transfer  rate  at  this  Interface  has  been  determined  from  an  empirical 
correlation  which  will  be  discussed  later. 

In  review,  the  basic  ablation  process  occurs  as  follows: 

When  the  wall  temperature  reaches  a  certain  level,  which  depends  on  the  surface 
material's  reactivity,  a  heterogeneous  chemical  reaction  (oxidation)  between  the  hot 
boundary  layer  and  the  surface  material  ensues.  In  this  first  "reaction-rate 
controlled"  regime,  the  surface  oxidation  rate  increases  exponentially  with  surface 
temperature  (e.g.,  by  an  Arrhenius  kinetic  mechanism).  Reaction-rate  control 
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persists  with  increasing  surface  temperature  until  the  diffusional  transport  of 
oxident  to  the  surface  becomes  too  slow  to  maintain  the  exponential  increase  in 
reaction  rate.  This  "transition  regime"  extends  to  the  temperature  level  where 
the  surface  reaction  becomes  completely  diffusion  controlled.  At  this  condition, 
the  reaction  is  said  to  be  in  the  diffusion-limited  regime  and  is  characterized 
by  surface  oxidation  occurring  in  stoichiometric  proportions  (the  stoichiometry 
depending  on  the  products  formed).  As  the  surface  temperature  continues  to  in¬ 
crease,  the  oxidation  rate  remains  diffusion-limited  until  a  critical  temperature/ 
pressure  condition  is  reached  at  which  the  graphitic  material  also  begins  to  sub¬ 
lime.  In  this  sublimation  regime  the  mass  loss  rate  again  increases  exponentially 
with  surface  temperature  as  dictated  by  the  partial  pressure  of  the  vaporizing 
carbon  species. 

For  our  combustor  ablation  model,  the  sublimation  regime  has  been  neglected'. 
This  simplification  has  been  assumed  because  sublimation  of  graphite  materials 
have  been  observed  to  begin  at  approximately  5500°R,  which  is  a  temperature  well 
above  the  maximum  expected  surface  temperature  under  consideration  for  the  tra¬ 
jectories  presently  being  evaluated. 


THE  REACTION-RATE  CONTROLLED  (KINETIC)  REGIME 

Experimental  surface  reaction  rate  data  are  normally  presented  in  the  Arrhenius 
plane,  i.e.,  log  m^  vs  1/T,  in  an  attempt  to  correlate  the  data  with  the  empirical 
Arrhenius  rate  equation. 


“r  “  (-E/RT^) 


(1) 


where  mj.  is  the  carbon  mass  loss  rate,  kg  Is  the  "frequency  factor"  or  "reaction 
rate  constant",  n  is  the  order  of  the  reaction,  E  is  the  activation  energy,  and 
is  the  partial  pressure  of  the  oxidant  at  the  surface. 

Scala^  in  an  attempt  to  bracket  the  available  data  on  graphite,  utilized  a 
set  of  "slow"  and  "fast"  values  for  E  and  k©  with  n  =  1/2.  Subsequently  Metzger^, 
based  on  the  experimental  arc- jet  data  on  ATJ,  formulated  a  third  set  that  resulted 
in  rates  between  the  "fast"  and  "  slow"  values  of  Scala  and  labeled  as  the  'moderate' 
rates.  The  resulting  values  are: 

k^  =  9.65  X  10^  lb/(ft^-  sec-atm^^^) 

.  E  =  44,000  cal/mol 


R  =  1.98718  ^Ral/mole  -  °K  (universal  gas  constant) 

T^  *  wall  te^erature,  °K 


Presently,  there  are  no  reported  experimental  data  on  carbon-carbon  to  arrive 
at  similar  reaction  constant  values  for  this  material.  Therefore,  the  above 
values  for  ATJ  graphite  are  employed  as  the  first  approximation  for  the  E  and  k© 
of  carbon-carbon  for  our  present  model. 

*  t...  \  fi  . 

Nagle  aiid  Strlckland-Constable  have  reported  rate  constants  for  the  oxida¬ 
tion  of  pyrolytic  graphite  (PG) .  They  fit  their  experimental  data  to  a  surface 
reaction  mechanism  proposed  by  Blyholder,  Blnford,  and  Eyrlng^.  The  mechanism  was 
developed  to  explain  the  observation  that  the  oxidation  of  PG  appears  to  follow  at 
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least  two  different  rate  equations.  At  temperatures  below  1000°K  (1340°F) ,  the 
oxidation  can  be  described  by  a  rate  equation  with  an  activation  energy  of  25 
to  30  kcal/gmole.  The  rate  of  oxidation  reaches  a  maximum  in  the  region  between 
1000  and  1700  K  (13A0  to  2600°F),  then  decreases,  and  at  still  higher  temperatures 
again  Increases  with  an  activation  energy  from  70  to  90  kcal/gmole. 

It  is  assumed  that  two  types  of  surface  active  sites  are  present  on  PG.  Type 
A  sites,  the  more  reactive  sites,  support  the  oxidation  reaction  at  low  temperatures. 
At  higher  temperatures  some  Type  A  sites  are  converted  to  less  reactive  type  B 
sites,  accounting  for  the  temporary  decrease  in  the  oxidation  rate  in  the  inter¬ 
mediate  range  of  temperatures.  This  mechanism  is  described  by  a  set  of  three 
elementary  reactions  wherein  CO  is  assumed  to  be  the  only  product: 


(1)  A  +  O2 

A  +  2CO 

(i) 

(ii)  B  +  O2  -»■ 

A  +  2C0 

m^  (ii) 

=  kg  (1-X) 

(iii)  A  B 

mj.  (iii) 

=  k^X 

where 

X  =  Fraction  of  surface  covered  with  Type  A  sites 

(1-X)  =  Fraction  of  surface  covered  with  Type  B  sites 

=  Partial  pressure  of  oxygen,  atm  (=  P^  * 

$^2  “  Mole  fraction  of  oxygen 

^A’^B’^z’^T*  constants 


\ 


where 


w 


The  total  rate  of  reaction  is 


“r 


X  +  kg  Pq 


(1-X) 


(2) 


(3) 


2 

in  gram  atoms  carbon/cm  -sec.  Values  for  the  reaction  rate  constants  were  obtained 
by  finding  those  values  which  best  fit  the  experimental  data: 

'■'iT 

k^  ■  20  exp  (-30,000/RT)  gram-atoms/ (cm  -sec-Atm)';-  . 


*^B 

kip 


2 

4.46  X  10“^  exp  (-15,200/RT)  gram-atom/(cm  -sed-AtmX-. 
1.51  X  10^  exp  (-97, 000/RT)  gram-atom/ (cm  -sec)  \  ,  ' 

21.3  exp  (+4,100/RT)  1/atm.  ,  , 
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TRANSITION  AND  DIFFUSION  LIMITED  REGIMES 


In  the  transition  regime  the  mass  loss  rate  is  computed  from, 


1 

2 


which  is  iwMMK  1/2-order  reaction  (Reference  1)  and  where  m^  is  the  ma 
rate  in  the  diffusion- limited  regime. 


The  evaluation  of  graphite  surface  ablation  throughout  all  the  thermo¬ 
chemical  regimes  is  summarized  by  the  non-dimensional  oxidation  term  called  the 
"mass  transfer  parameter" ,S,  which  is  defined  as 

3  =  m/puC^ 

where 


m  =  the  mass  loss  rate 

puC^j  =  heat  transfer  coefficient 

Cj^  =  Stanton  N^imber 

and  for  the  diffusion  limited  regime 

6p  ==  mj)  /(puCj^)  (6) 

Equation  (6)  can  be  re-arranged  to  obtain  the  flip  from  the  "diffusion  limited 
mass  transfer  parameter".  Bp 


ip  *  (puCj^)  6^  (7) 

where  (puCj^)  and  Bp  can  be  Independently  determined,  which  will  be  discussed 
later. 

The  mass  loss  ra^-e  equations  for  each  of  the  thermochemical  regimes  are 
combined  by  means  of 'the  Mickley-Spaldlng  blowing  relation®  to  form  the  following 
equation  to  determine  B’  (which  is  essentially  the  solid  carbon  to  free  stream  gas 
mixture  mass  ratio) , 


(8) 


REACTIVITY  OF  CARBON 

At  elevated  temperatures,  and,  to  a  first  ajpproxlmatlon,  irrespective  of 
the  nature  of  the  carbon,  the  reactivity  is  the  greatest  in  oxygen  and  least  in 
hydrogen,  as  summarized  in  Table  1  (Reference  9).  This  table  also  reveals  that 
the  C  +  62  reaction  is  roughly  10^  times  as  fast  as  the  C  +  CO2  reaction  at  around 
1900°R.  Based  on  this  finding,  our  model  assumes  th(|Lt  the  carbon  wall  reacts  only 
with  the  O2  remaining  in  the  by-products  of  the  Shell^yne-H/alr*  mixture  and  CO 

*Shelldyne-H  is  the  fuel  under  consideration.  See  Appendix  A  for  properties. 


is  the  only  by-product  of  the  solid  carbon  reaction  according  to 


2C  +  ^  2C0 

Table  1 


Approximate  Relative  Rates  of  Carbon  +  Gas 
Reactions  at  193l‘^R  and  0.1  Atm 


Reaction 


Relative  Rates 


(9) 


C  +  CO2 
C  +  O2 
C  +  H2O 
C  +  H2 


1 

1  X  10^ 

3 

3  X  10"^ 


DIFFUSION-LIMITED  MASS  TRANSFER  PARAMETER  (B^) 

The  diffusion-limited  mass  transfer  parameter,  Bq,  is  essentially  the  solid 
carbon  to  free  stream  gas  mixture  mass  ratio  at  stoichiometric  proportions.  With 
the  assumption  that  the  solid  carbon  only  reacts  with  the  O2  in  the  free  stream 
and  CO  is  the  only  by-product,  the  basic  reaction  equation  is 


C  + 


2(0.21  -  0.198ER) 


“a 

^  ‘fg'ER'Cj^^  H^g  +  0.198’ER*02  +  0.79*N2  + 


(0.21  -  0.198*ER)  O2  -*■  CO  +  (Shelldyne  H/alr  less  O2)  (10) 

The  relationship  defining  Bq  as  a  function  of  the  fuel/air  equivalence  ratio 
(ER)  can  be  derived  with  the  utilization  of  this  reaction  equation.  After  further 
rearrangement  of  terms  and  the  substitution  of  molecular  weights,  the  Bq  equation 
reduces  to 


,  (0.17A8  -  0.1648  ER) 

Bd  •  (1*0  +  0.0730  ER) 


(11) 


Figure  1  presents  the  variation  of  Bd  versus  ER  per  equation  (11) .  When  the  ER 
equals  zero,  the  Bj)  is  essentially  equal  to  the  diffusion  limited  value  for  air 
(0.1748). 


COMBUSTOR  HEAT  TRANSFER  CORRELATION 

Over  many  years  of  combustor  testing  at  APL,  heat  transfer  measurements  have 
been  taken  for  varied  geometries  and  different  fuels.  It  has  been  the  intent  of 
this  data  gathering  exercise  to  correlate  the  average  heat  flux  values  from  these 
tests  against  a  heat  release  parameter  which  is  proportional  to  the  effective 
equivalence  ratio  (ER^ff)  for  any  given  fuel.  This  heat  release  parameter  is 
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equal  to  the  product  of  the  fuel-air  ratio  (f),  the  lower  heating  value  of  the 
fuel  (hf),  and  the  combustor  efficiency  (n) . 


Figure  2  (from  Reference  10)  presents  the  correlation  of  the  normalized 
wall  heat  fliuc  parameter  as  a  function  of  the  heat  release  parameter  (the  ER^ff 
scale  for  Shelldyne-H  is  also  Included.)  The  normalized  wall  heat  flux  para- 


meter  is  equal 

to 

/  /.  V  /  w  Ah 

/  i  — § — 

■  A 

; 

where: 

Wall  heat  flux 

A 

Wall  area 

V  .  ..  ^ 

Aci 

Combustor  inlet  area 

Wa 

= 

Air  mass  flow  rate 

Ah 

- 

Average  gas  to  wall  enthalpy  difference 

Notlng  that  by  definition 

PuCj^  »  (q„/A„)/Ah  (12) 

the  combustor  heat  transfer  coefficient  can  be  determined  for  any  air  mass  flow 
rate  (Wg)  and  a  specific  combustor  geometry  cross  sectional  area  (Aj.3  from  this 

^  11 

empirical  correlation  instead  of  an  analytical  approximation  such  as  Bartz 
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ERgpp,  Shelldyne-H 
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<  0.0020 


S  0.0008 


6  8  10  12  14  16 

Heat  release  parameter  (f  •  •  hf)  x  10"2  BTU/lb 

Fig.  2  Correlations  of  wall  heat-flux  parameter  with  heat-release  parameter. 
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THERMOCHEMISTRY  DATA  DETERMINATION 


The  thermochemlcal  erosion  of  graphite  materials  Is  dependent  on  Information 
denoted  as  thermochemistry  data.  This  data  consists  of  the  mass  transfer 
parameter  (3')  and  the  enthalpy  of  the  gas  mixture  at  the  graphite  wall  Interface 
(hw)  tabulated  for  different  wall  temperatures  (T^)  and  pressures  (?«) •  The 
ablation  model  reviewed  In  the  previous  sections  Is  now  utilized  to  calculate 
the  3'  for  varying  T„  and  ^w  for  a  specific  supersonic  combustor  design  and  fuel 
equivalence  ratio  (ER) . 

A  basic  procedure  leading  to  the  calculated  3'  is; 

1.  For  the  desired  ER,  the  3^  and  the  oxygen  mole  fraction  of  the  by¬ 
products  mixture  can  be  determined  for  the  Shelldyne-H/air  fuel 
mixture  (See  Appendix  A  for  Shelldyne-H  properties) ; 

2.  For  the  specific  combustor  design  (A,  M,  Tg) ,  the  air  mass  flow  rate 
can  be  determined  for  varying  Pw* 

3.  The  heat  transfer  coefficient  between  the  free  stream  gas  mixture 
and  the  combustor  wall  can  then  be  determined  from  the  heat  flux 
correlation  for  the  ER  and  w^; 

4.  The  mass  loss  rates  for  each  of  the  thermochemlcal  regimes  (depending 
on  PG  or  C/C)  is  determined  for  varying  T^,  and  ^w» 

5.  3*  can  be  determined. 

The  relationship  between  3 ' ,  h^,  Tw»  is  generated  by  means  of  a  thermo¬ 
chemlcal  computer  code  (NOTS)^^.  The  calculations  are  based  on  equilibrium 
thermochemistry  where  solid  carbon  is  the  fuel  reacting  with  the  hydrogen, 
carbon,  nitrogen,  and  oxygen  constituents  of  the  gas  mixture.  Twenty-one 
species  have  been  considered  In  the  equilibrium  calculation  In  which  their 
thermochemlcal  data  are  from  the  JANNAF  tabulations  where  all  enthalpies  are 
referred  to  the  base  at  298.16°K,  1  atm  and  the  elements  In  their  natural  state. 


COMPUTER  PROGRAM 

The  above  procedure  has  been  Incorporated  Into  a  computer  program.  A  new 
version  of  the  NOTS  (N0TS-CMA)12  code  has  been  created  with  modifications  and 
the  addition  of  a  subroutine  for  the  sole  purpose  of  generating  the  thermochemistry 
data  In  the  proper  Input  format  for  the  heat  transfer  analysis  computer  program 
CMa13, 


Figure  3  presents  the  calculated  3'  for  various  Shelldyne-H  ER's  over 
carbon-carbon  material.  In  each  case,  the  3'  levels  off  at  the  diffusion 
limited  value.  For  any  specific  ER,  the  effect  of  P^  on  3'  is  negligible 
above  2500°R.  However,  for  the  case  of  Shelldyne-H  over  PG,  its  an  entirely 
different  result.  Figure  4  presents  the  3*  for  the  ER  0.5  case  showing,  not 
only  a  variation  with  wall  pressure,  but,  more  Importantly,  the  mass  transfer 
parameter  reaches  a  maximum  In  the  temperature  range  of  3500-4000°R  then  drops. 
The  oxidation  of  PG  Is  klnetlcally  controlled  (not  diffusion  limited)  and  behaves 
In  a  peculiar  manner  per  the  previous  discussion.  Since  CMA  requires  that  the 
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6'  versus  T„  be  monatonically  Increasing,  6'  is  bounded  at  a  maximum  and  min¬ 
imum  slope  as  shown  in  Figure  4  and  a  range  of  surface  ablation  is  calculated. 


EXPERIMENTAL  RESULTS 

The  development  of  high  performance  supersonic  combustor  materials  re¬ 
quires  an  understanding  of  the  mechanical  and  chemical  processes  taking  place 
at  the  ablating  surface  of  various  materials.  This  is  impossible  to  obtain 
from  analysis  alone  and  requires  some  experimentation.  Consequently  a  series 
of  tests  were  performed  at  the  McDonnell  Douglas  High  Impact  Pressure  (HIP)^^ 
arc  heater  facility  to  simulate  sea  level  and  high  altitude  cruise  conditions 
in  a  supersonic  combustor. 

The  purposes  of  these  tests  were  to  evaluate  candidate  supersonic  com¬ 
bustor  materials  and  to  obtain  mass  loss  and  surface  temperature  data  for 
comparison  with  analysis.  A  total  of  40  models  were  tested,  with  each  model 
nominally  0.5  inches  wide  by  0.5  Inches  thick  by  1.25  inches  long.  A  model 
was  placed  on  each  arm  of  an  eight  arm  turret  (Figures  5  and  6)  and  injected 
into  the  arc  jet.  Thirty-two  of  the  models  were  tested  at  a  10°  angle-of- 
attack  for  75  s  each  at  a  nominal  brightness  temperature  of  5000°R  (e=l) 
using  a  similated  Shelldyne  -  air  mixture  to  approximate  high  altitude  cruise 
conditions.  Eight  models  were  tested  for  40  s  at  a  25°  angle-of-attack  using 
a  Shelldyne  -  air  mixture  to  simulate  sea-level  conditions.  Surface  brightness 
temperature  was  measured  at  four  locations  along  the  model  centerline  (0.25, 

0.55,  0.85  and  1.15  inches  from  the  leading  edge)  using  optical  pyrometers. 

The  model  backface  temperature  history  was  measured  at  the  0.85  inch  position 
using  a  tungsten-rhenium  thermocouple.  Model  recession  histories  were  obtained 
from  16  mm  movie  films  of  the  model  profile. 

Several  materials  were  used  to  determine  which  were  best  suited  for  scrara- 
jet  engine  combustors.  Particular  Interest  was  given  to  metal  impregnated  carbon- 
carbon  (C-C)  materials  which  in  theory  form  protective  metal  oxide  layers,  re¬ 
tained  at  the  surface  by  high  viscosity  silica  compounds  also  added  to  the 
specimens.  A  summary  of  the  processing  treatments  and  composition  of  metal  and 
silica  additives  for  each  specimen  is  presented  in  Table  2.  In  addition,  some 
reference  materials  Including  virgin  C-C,  pyrolytic  graphite  (PG)  coated  C-C, 
pyrolytic  graphlte/silicon  carbide  (PG/S-C)  coated  ATJ  graphite,  and  ATJ-S 
graphite  were  also  used. 

A  summary  of  the  mass  recession  loss  rate  results  is  presented  in  Figure  7 
for  the  long-range,  high  altitude  tests  and  Figure  8  for  the  sea-level  test.  In 
the  high-altitude  tests  there  is  no  particular  trend  in  the  metal  additive  coirpo- 
sltlon  or  process  and  there  appears  to  be  no  benefit  from  metal  additives.  On 
the  contrary,  the  limited  data  suggests  metal  additives  have  increased  the  re¬ 
cession  rate.  Except  for  model  #21  (ZR(.6)  -  Si(.2)  -  B(.2)/RH),  the  pyrolytic 
graphite  is  clearly  superior. 

The  recession  rates  for  materials  tested  at  the  sea-level  condition  show 
that  the  metallzed  C-C  samples  (24,  14  and  2)  are  the  same  or  worse  than  the 
plain  C-C  sample  number  45.  Moreover,  PG  and  PG/SiC  had  much  lower  recession 
rates.  Although  it  appears  that  the  metallzed  C-C  materials  offer  no  significant- 
improvement  over  other  combustor  materials,  shear  stresses  on  the  surface  were 
slightly  higher  than  desired  and  may  have  altered  the  performance  that  these 
materials  would  experience  in  a  real  combustor  application. 


m 


m 


Table  2  MaterTkl  Specimen  Number  Identifications  for  HIP  Tests 


4-  4-4* 

Prora 

ssing  Precursor  an 

d  Phase 

vi  1  lie  Additives  Composition 

RH 

lit 

Jet 

AO 

Zr(.f))-ai(.2)-B(.2) 

15 

s> 

<g) 

17 

22 

iE 

38 

40 

Zr(.8)-Si(.l)-B(.l) 

0 

1 

14 

1 

1 

Zr(.6)-Si(.2)-Ti(.2) 

Zr(.8)-Si(.l)-Ti(.l) 

) 

Hf (,6)-Si(.2)-Ti(.2) 


N'oa-netallized  specimens  (pitch  only) 


are  @  ^  ^ 


Cijde 


Vi  -•  .iLkoxy  -  HiC 

-  A i  ■■•cozy  -  LoC 

■<K  -  Ras inace  -  HiC 
>  =  Re^-inate  -  I,oC 

-  Reninate  -  Oxide 
‘•0  =  Alkr.xy  -  Oxide 


Runs  1554  &  1555 
Runs  1558  &  1559 
Run  1564 


'^"'•'iribars  'n  parentheses  refer  to  tr.ole  fraction  of  the  metal  additive. 

! nl.s  ire  listed  from  top  to  bottom  in  order  of  incrocising  theoretical 
•'  r  n. ;  oriness . 
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21  PG 
Fig.  7 


37  13  47,48  23  16  9  1  ATJ-S  7  39 

Material  Designation 

Material  recession  rate  performance  for  high-altitude  tests. 


18  ,  '25 
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PG/SiC  PG  Graphnol  45  ATJ-S  24  14  2 

Material  designation  ' 

Fig.  8  Material  recession  rata  performance  for  sea-level  tests. 
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COMPARISON  OF  MEASURED  AND  PREDICTED  MATERIAL  RECESSION 

rLScal  properties  of  carbon-carbon  were  obtained  from  the  EST  computer  codeia 
assuming  no  CO2  in  the  chemically  reacting  products  near  f 2C0) 
ThJ  rational  for  this  assumption  is  based  on  the  fact  that  the  (2C  +  0^2 
faction  ;fte  is  much  fasted  than  the  (2C0  H-  O2  -  2C0) 

turn,  is  much  faster  than  the  (2C0  +  O2-  2C0)  rate.  ’ /“haLf 

duced  and  then  swept  away  from  the  reacting  surface  before  it  has  a  chance 

to  react  with  02»  to  form  C02* 


The  analytical  model  consists  of  a  0.5  Inch  thick  piece  of  ATJ-S,  insulated 
from  a  0.5"  steel  support  by  a  0.02"  layer  of  zirconium  dioxide  (Figure  9). 
Measured  surface  temperatures  were  obtained  from  black  body  optical  pyrometer 
temperatures  assuming  an  ATJ-S  surface  emissivity  of  0.78  at  0.9p. 

Predicted  surface  temperatures  compare  favorably  with  optically  measured 
temperatures  during  the  first  20  to  30  seconds  of  the  test  (Figure  10).  After 
this  time  however,  measured  temperatures  decrease  sharply  and  are  no  longer  in 
good  agreement  with  predicted  temperatures  which  continue  to  Increase.  The 
reason  for  this  is  that  as  the  specimen  ablates  it  recedes  out  of  the  freejet 
main  stream  and  into  the  cooler  bovindary  layer.  Therefore,  only  the  first  20 
to  30  seconds  of  the  test  are  a  valid  simulation  of  actual  combustor  conditions, 
and  data  taken  after  that  should  be  ignored. 

Surface  ablation  rates  measured  during  this  test,  compare  well  with  predicted 
recession  rates  throughout  the  test  (Figure  15 .  Both  experiment  and  analysis 
Indicate  ablation  begins  between  6  and  11  seconds  and  continues  at  a  constant 
rate  throughout  the  test.  Constant  recession  rates  of  0.0027  and  0.0024  inches/ 
sec  have  been  measured  and  compare  with  a  slightly  higher  predicted  rate  of  0.003.< 
Inches/sec.  The  somewhat  surprising  aspect  of  this  data  is  that  measured  re¬ 
cession  rates  remain  constant  even  at  times  greater  than  30  seconds,  when  the 
surface  temperature  is  decreasing  and  the  heated  surface  is  moving  into  the  free¬ 
jet  boundary  layer.  This  seemingly  Inconsistent  behavior  can  be  explained  as 
follows:  In  general,  ablation  rate  is  expressed  as  8'  *  PuCi^  (Equation  7), 
where  8'  is  the  diffusion  mass  flux  ratio  expressed  as  a  function  of  the  material, 
its  surface  temperature  and  pressure,  while  puC^  is  the  enthalpy  based  heat 
transfer  coefficient.  For  ATJ-S  the  ablation  rate  is  diffusion  limited  at 
temperatures  above  2400°R,  which  means  8*  is  constant  about  2400°R.  Since 
surface  temperatures  remain  above  2400°R  for  most  of  the  test  (Figure  10),  8' 
is  a  constant.  However,  in  order  to  have  a  constant  recession  rate,  puC]^  must 
also  be  constant.  Apparently,  as  the  surface  recedes  into  the  cool  thermal 
boundary  layer,  the  effect  of  decreasing  stream  velocity  (u)  is  offset  by  a 
decrease  in  the  reference  length,  since  the  effective  distance  to  the  leading 
edge  is  decreasing  as  the  leading  edge  leaves  the  freejet  main  stream.  The 
good  agreement  between  experimental  and  analytical  surface  temperatures  and 
ablation  rates  give  us  added  confidence  in  our  analytical  model. 
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Qaero 


Zirconium 

ATJ-S 
(25  nodes) 

dioxide  ^ 

(1  node) 

Steel 

(10  nodes) 

Analytical  model 


Fig.  9  HIP  test  specimen. 
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Appendix  A 

Shelldyne-H  Fuel  Properties 


The  raajet  fuel  under  consideration  for  the  hypersonic  missile 
Is  the  hydrocarbon  fuel,  Shelldyne»H  (Cx4  ^18^*  of  Its  basic 

properties  are: 

Density  ■  67.5  #/ft^ 

Stoichiometric  fuel-air  ratio  (fg)  >  0.072832 

Lower  heating  value  ■  17,890  Btu/# 

Table  A-1  presents  a  tabulation  of  the  mole  percent  of  various  by¬ 
products  of  combustion  for  varylxig  equivalence  ratio  (ER)  of  the 
Shelldyne-H/alr  mixture.  The  mole  percent  of  O2  by-product  versus 
ER  Is  plotted  for  two  different  free  stream  air  total  temperatures  of 
1080*R  and  3780*R  In  Figure  A-1.  The  approximate  curve  fit  of  oxygen 
mole  function  for  use  In  our  analytical  model  Is 

(402)moles  ■  0.21  -  0.198  ER 
which  Is  also  plotted  In  Figure  A-1  for  comparison. 

Figure  A-2  presents  the  variation  of  combustor  total  gas  temperature 
versus  free  stream  air  total  teoqierature  for  varying  equivalence  ratio. 


11-124 


Combustor  total  temperature  (‘’R) 
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Free  stream  total  temperature  (°R) 

Fig.  A-2  Combustor  versus  free  stream  T^,  Shelldyne  •  H. 
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